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Abstract 



In the framework of the conjectured duahty relation between large gauge theory 
O I and supergravity the spectra of masses in large A^ gauge theory can be determined 

^ ' by solving certain eigenvalue problems in supergravity. In this paper we study the 

eigenmass problem given by Witten as a possible approximation for masses in QCD 
D , without supersymmetry. We place a particular emphasis on the treatment of the hori- 



zon and related boundary conditions. We construct exact expressions for the analytic 



^ ■ expansions of the wave functions both at the horizon and at infinity and show that re- 

quiring smoothness at the horizon and normalizability gives a well defined eigenvalue 
problem. We show for example that there are no smooth solutions with vanishing 
derivative at the horizon. The mass eigenvalues up to = 1000 corresponding to 
smooth normalizable wave functions are presented. We comment on the relation of 
our work with the results found in a recent paper by Csaki et al., [hep-th/ 9806021 
which addresses the same problem. 
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1 Introduction 



The problem of solving QCD in the nonperturbative large limit has been out- 
standing for several decades [|l], §]• It has been suspected that the answer to this 
question will come from string theory. Recently a very interesting proposal 0] has 
been introduced and further explored |§, 0, § which involves a relationship between 
large N super Yang-Mills theory and AdS supergravity |^-|2^. This correspondence 



which was first investigated in studies of 3-branes gives the possibility of studying 
large N properties of Yang-Mills theories using classical supergravity 0. The later 
is expected to give results that should be valid for the strongly coupled gauge theory. 
At present time comparison of the two theories was done for operators and correlators 
protected by supersymmetry [|T^-[^. Other quantities like the entropy or Wilson 
loops represent predictions of the conjecture 0]-p5i- general systems involving 



p-branes a notion of generalized conformal symmetry was found in ||2^. One can 
expect that a similar correspondence holds also in theories without supersymmetry 
and ultimately in QCD. Witten has presented such an extension where proprties 
of finite temperature Yang-Mills theories are to be computed using AdS black hole 
backgrounds in gravity 0. 

According to Witten's generalization of the conjecture by Maldacena in in 
order to study A/" = 4 super Yang-Mills theory at large A^, high temperature and 
strong t'Hooft coupling, one should consider the Euclidean Schwarzchild black hole 
solution in AdS^ x space-time in the limit where the black hole mass is large [Q. 
In this limit the metric can be written as 

dr'^ 



\ f)2 j,2 J 1=1 



where dflf is the round metric on 5*^, r = 6 is the horizon radius and the coordinate 
r is the Euclideanized periodic time coordinate. This metric is obtained as a solution 
to the type IIB supergravity equations of motion following from the 7 — »• limit of 



the action 

S = / d'^x^iR - + ■ ■ ■ + jeM-lm + ■■■) (2) 

where the order 7 = l/8({3)a'^ terms contain the first string corrections to super- 
gravity and where is a certain combination of terms quartic in the Weyl tensor 
P^ , The inclusion of stringy a corrections corresponds to including strong cou- 
phng expansion corrections in the gauge theory ^ |[. 

In this paper we would like to study the proposal in ^ for the supergravity 
calculation of the mass gap in QCD. In the next section we will examine the equations 
of motion for free scalar field propagation on AdS black hole backgrounds by rewriting 
them in the form of an Hamiltonian problem. We will then address the problem of the 
behaviour of the wave function at the horizon. In section 3 we present our exact results 
for the wave functions and show that there are no normalizable smooth solutions with 
vanishing derivative at the horizon. Using the exact form of the solutions we then 
exhibit the glueball mass eigenvalues predicted by supergravity. Finally, in section 4 
we close with some conclusions. 



2 Black Hole Backgrounds 



According to Witten the equations for free field propagation, d^^^g^^dyr]) = for a 
scalar field, in the five-dimensional space-time described by the first terms of (0) with 
r compactified on should give glueball masses for QCD^. Similarly a computation 
in an AdSj black hole background is expected to be of relevance for glueball masses 
in QCD4 p|. One should look for solutions behaving like (static) plane waves along 
the Xi directions ~ ri{r)e^^'^ and then demanding normalizability and regularity of 
the behaviour of ri{r) at r = b and r = 00 will select only certain allowed discrete 
values for = —k"^. These values of are then interpreted as particle masses 



3 



in the three-dimensional world parametrized by the Xj. To study the corrections to 
these masses in the strong coupling expansion one then should work with the 0(7) 
corrections to the background metric (|lD and to the dilaton field. To 0(7) and for 
the purpose of computing mass corrections, it is consistent to take a classical solution 
with a vanishing dilaton field. The 0(7) correction to the metric was found in |^| 
and one uses it to compute the 0(7) corrections to the glueball masses. 

Consider the general metric for the AdSn+i black hole in the large mass limit 

i 

2 Ln-2 J 2 n-1 

" ' 1^2 i=l 

The equation of motion for a free scalar field of the form ~ jj^j-^j^^k-x 
given by 

drir'^''^-^ - ^)9rV) + r^-'m'v = (4) 

where = —k"^ is the {n — l)-dimensional mass. Consider the measure coming from 
the metric (^) above (we set 6 = 1 for the remainder of this section) 

< r]\r] >= / drr"~^r](r)r]* (r) . (5) 



1 

In order to trivialize the measure we can take a new variable y = for which with 
^{y) = y^viy) the equation becomes 

d,{{y' - l)dy<^) + (--(3 + y-') + -^y-")<^ = (6) 

Integrating from the horizon to 00 against and integrating by parts assuming 

normalizability and smoothness we obtain a bound on the possible values of m^. For 
example, for = 4 we find that > 4. To eliminate the first derivative term 
in and to write that equation in terms of an Hamiltonian problem we now take 

TL 1 _ 

y = r2 = cosh(w) and redefine A{w) = smh{2w) 2 ri{w) . This gives 

IdlAiw) - Viw)Aiw) = (7) 



where the potential is now given by 
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2V{w) = 1 — sinh(2tt;)~^ r-m^ cosh(w) 



We are interested in the wave function A{w) for the zero eigenvalue of (0). If we 
expand the potential around the horizon w = we obtain V{w) = — 1/ {Sw"^) + (2/3 — 
{2m?)/n^) + ((4m^)/'n,^ — 2/lb)w'^ + O(w^). The harmonic oscillator perturbed by 
a potential of the form A(l/w^) was examined in |26]. Our potential corresponds 
precisely to the limiting case A = —1/8 in that reference beyond which the Hamil- 
tonian is not bounded below. The indicial equation for (|^ with the potential ex- 
panded about the horizon w = 0, will have a double root \. Therefore, near the 
horizon we will have the behaviours A{w) ~ and A{w) ~ log(w) for the two 
independent solutions of (|^). Both solutions are normalizable near w = Q and we 
also have normalizable density of probability currents at the horizon of the form 
J{w) ~ A{w)du,A{w) ~ const, or ~ log^(w). The two solutions give wave functions 
?7(r) for which behave near the horizon like rjiw) ~ const, or rjiw) ~ log(w). The 
first derivatives then become drj/dr ~ const, or drj/dr ~ const. /w. We therefore 
expect that the Neumann boundary condition may never be attained at the horizon 
for a regular solution. Indeed we note that we have a potential which is singular at 
the horizon and that it could be expected that it is not possible to demand Dirichlet 
or Neumann boundary conditions there and as we will see this is what happens in our 
case. It would be interesting from this general point of view to understand if possible 
tunneling effects could contribute in a small amount to the values of the eigenmasses. 
Our solutions of equation (|) which we will present in the next section are consistent 
with the above behaviour. 

To formulate the eigenvalue problem, one fixes the behaviour at oo such that the 
solution is normalizable. Then demanding regularity of the solution at the horizon 
determines a discrete set of masses. The equations that describe the wavefunctions 
corresponding to motion in the AdS black hole backgrounds have regular singular 
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points at 0, 1, horizon and oo and also at other points according to the value of n. 
In view of the discussion above, it might be tempting to ask for solutions which are 
regular at the origin (instead of the horizon) and which decay well enough at oo, 
and hope that this would define an interesting eigenmass problem. However, a closer 
look at shows that the eigenvalues even if they exist are not guaranteed to be 
positive in that situation. 



3 Calculations and Results 

To leading order in 7 the equation of motion for the quadratic fluctuation 7]q of the 
dilaton field is § 

9r.(r(r^ — b'^)drrio) + mlb'^rrjo = (9) 

where one takes b < r < 00 and where mo is the leading contribution to the mass 
in the strong coupling expansion. The eigenvalues ml will provide the masses of the 
scalar glueball states. Considering first the behaviour of the solution at infinity 
it is useful to rewrite the equation in the variable z = b/r with < 2 < 1, 

One wants to find normalizable wave function solutions of @ and this fixes the 
behaviour at 00 to be like tjq ~ This behaviour at 00 provides us with 

a Taylor expansion for rjo around 2; = of the form 770 = I]^2 Cn-z^" where to fix 
the overall normalization of rjo we take ci = 0, C2 = 1 and then obtain the recursion 
relation for n > 2 

Cnjml/b^) - c„_i(2(n - l)(2n - 3) + 2(n - 1)) 

(2n + 2)(2n+l)-6(n + l) ' ^ ' 



We next concentrate on the behaviour of the solutions of this equation near the 
singularity at the horizon. We will first find an expression for the analytic solution 
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at the horizon and use it to reduce the order of the equation and show that the other 
independent solution is not smooth at the horizon. In order to better describe the 
vicinity of the horizon let us use the variable ( = b'^/r'^ — \ = \jz^ — \ such that the 
horizon is at = 00. The equation becomes 

(i^rin ,1 1 1 , drin rrin , 1 1 1 < , ^ 

where we can expand the fractional coefficients in powers of C and where we take a 
power series ansatz ?7o(C) = S^o^nC"- We obtain that the first coefficient 6o is free, 
hi = (mQ/86^)6o and that the other coefficients can be determined in terms of Bq from 
the recursion relation 

k — fc — 

(13) 

At the horizon this solution goes to a constant bo and the first derivative drjo/dr = 
{-2/b)dr]o/d( = {-2/h)hi = {-ml/4:b^)bo ^ 0. We will now use this solution to 
reduce the order of the equation and find a second linearly independent solution. Let 
be the analytic solution defined by the recursion relation ( |T3| ) and set tiq = ip ■ g. 
Inserting in ([T^) we obtain a first order equation for the ^-derivative of g, 

d'^q,dq^ -I dih ^ , 1 1 1, , , 

implying that dg/dC, = const. ((^ + !)/{({( + 2)^/'^). Therefore the second solution to 

(P^) has a first derivative which blows up at the horizon and is not smooth there. 

Consequently, there are no smooth solutions with vanishing first derivative at the 

horizon. This behaviour in the region close to the horizon is consistent with the results 

of section 2 where we have seen that the problem reduces to a Schrodinger problem 

for the harmonic oscillator perturbed by a potential of the form V{x) = (— f/8)l/x^. 

^Onc could take 1/r — 1 as well but it turns out that — 1 provides a much better behaviour of 
the coefficients of the power series for the regular solution and this is important to ensure a proper 
numerical treatment of the problem. 
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We would now like to use the exact form for the series solution at the horizon 
and try to fix the overall normalization by choosing the coefficient 60 in such a way 
that this solution matches with the Taylor expansion (|ll|) obtained by expanding at 
spatial infinity z = 0. Of course, such a matching of the two Taylor expansions over 
an interval will be possible only for certain values of rriQ and these are the values for 
the masses. We have used these exact expressions for the analytic expansions of the 
wave functions and have evaluated them numerically. We found that in practice this 
method yields strong conditions on the allowed values of itlq which can be found with 
arbitrarily high numerical precision. Indeed we found that for the allowed values of 
rriQ and once we compare the values of the two Taylor expansions at one point to fix 
the coefficient bo, the two Taylor expansions actually agree to a very high accuracy 
over an entire interval thus providing an impressive test of the method. Small changes 
in the values of mg away from the correct value are easily detected by the mismatch 
they produce between the Taylor expansions at the horizon and at 00. As an example 
we show the wave function for the 12-th eigenvalue = 895.8 in Fig.l below. We 
plotted the solution in x = b/r. The curve starting at the origin is determined from 
the Taylor expansion at infinity which we take up to a; = 0.9. We find that the power 
series converges extremely well in this region. From x = 1 we use the Taylor series 
from the horizon down to x = 0.8 where again we find a rapid convergence. As is 
clear from the figure the two expansions match perfectly as expected since we are 
describing the exact analytic form of the solution. We note that the radial derivative 
at the horizon is not zero (it is a factor of order one times the x-derivative) and is in 
fact not small. We find similar results for the other mass eigenvalues. 
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X 10' 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



Figl: The exact wave function for rriQ = 895.5. The horizontal axis is 
X = b/r. Prom a; = up to a; = 0.9 we use the Taylor expansion at oo. Prom 
X = 1 down to X = 0.8 we use the Taylor expansion about the horizon located 
at a; = 1. 



In Table I below we reproduce the first twelve values of ttIq. 
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SUGRA ml 
11.5877 
34.5270 
68.9750 
114.9104 
172.3312 
241.2366 
321.6265 
413.5009 
516.8597 
631.7028 
758.0302 
895.8410 



Table I 



The exact eigenvalue masses for the O^^ glueball in QCD^ derived 
from supergravity. Note that these exact supergravity masses have been 
rounded to the accuracy shown. 



The authors of used a "shooting" technique and numerically integrated the 
differential equation using the Taylor expansion at oo as an initial condition. To 
fix the values of ttIq one needs to fix the boundary condition at the horizon. If one 
uses the Neumann boundary condition r^Q = as proposed in p|, one finds numerical 
values for mp that are in excellent agreement with the values in Table I above even 
though there are no smooth solutions satisfying that boundary condition. Although 
the dependence of the eigenmass values on the boundary condition at the horizon is 
relatively weak, this is of course not true for the wave functions since one is precisely 
discussing the first derivative at the horizon. The fact that the actual eigenvalues 
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turn out to agree is interesting. It can be explained by the fact that in the "shooting" 
technique (&) is a rapidly varying function of mg as can be seen from Fig. 2 below. 
We can expect that the discrepancy bcwtcen the exact mass and the one obtained 
using the Neumann boundary condition will increase with increasing masses since the 
exact boundary condition has ?7q at the horizon increasing with ml. This point should 
be taken in consideration in future work on the subject and in particular in future 
studies of mass spectra from supergravity where one could conceivably demand a high 
accuracy. 



0.15 




_Q 25 I I I I ^ I I I I I I 

'11.58 11.582 11.584 11.586 11.588 11.59 11.592 11.594 11.596 11.598 11.6 



Fig2: This plot shows the dependence of the derivative at the horizon -^{b) 
as a function of mg in the "shooting" technique. The wave function is normal- 
ized so that the first term in the Taylor expansion about 00 is 1. In this 
normalization the exact wave function has ^ = —0.03 which would put the 
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exact mass at mn above 11.588 whilst the Neumann boundary condition would 



n2 



give ttiq below 11.588. 

We emphasize that our construction is based on matching the analytic forms of 
the wave functions over an extended interval. Consequently, the size of the error in 
the determination of the masses in our approach leads us to exclude the possibility 
of an exact mass formula of the form ml = 6n{n + 1). 

To find the 0(7) corrections to the masses of the 0^~^ glueballs one needs to 
study the equations of motion for quadratic fluctuations of the dilaton field in the 
metric background (|ip corrected to leading order in a . This correction was found in 
p5| where one can also find the expression for W. One sets = rriQ + 7771^ + 0(7^) 
and ri{r) = ?7o(r) + 7r7i(r) + 0(7^) and perturbs (H) about the vanishing dilaton 
background^. This gives the equation of motion (with 6=1) 

— (r(r - I))——) + rmQr]i = (15) 
ar dr 



(r 



. 300 600 1980, drio , . 2/75 75 165, 405 2I2O. 



Normalizability and the already known behaviour of ?7o once again fix the behaviour 
of ?7i at 00. By an analysis similar to the one we performed above one can show 
that there are no smooth solutions with vanishing derivative at the horizon. In this 
case the expressions for the Taylor expansions become a bit cumbersome. We have 
calculated the first few mass corrections vn\ by the method of matching the two Taylor 
expansions for the exact solutions. Once again for the same reasons that we explained 
above, we found values of vn\ which coincide with the ones obtained by Csaki et al. 
p8| via the "shooting" technique and we will not repeat those values here. Therefore, 
we also have nothing to add to the physical analysis that was done in that reference 
and in particular on the comparison with the results from the lattice [^7]]. 

■^Evcn though the dilaton is corrected to 0(7) as was calculated in ||2^ this does not affect the 
equation of motion for the fluctuations r\ to this order. 
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In |^8| the authors also studied the spectrum for the O glueball in three- 
dimensional QCD and the glueball mass spectra in four-dimensional QCD obtained 
from the black- hole geometry in AdSj x 5*^ P). The results of our analysis apply 
equally well to those cases. We note that in all these cases there is only one smooth 
solution at the horizon as the indicial equation always has a double root. This is 
physically interesting since otherwise matching with the behaviour at infinity would 
most probably not put enough restrictions on the wavefunctions and the eigenmass 
problem would likely be ill defined. 



In the case of the O glueball in QCD-^ one has the eigenvalue problem pS 



r(r^ - 1)^ + (3 + r^)^ + imlr - 16r=^)r/o = 0. (16) 

Normalizability fixes the behaviour at oo to be rjolr) ~ and the Taylor expansion 
at oo has the form ?7o(r) = I]J^Lo '^"^"^^ with x = 1/r^ (we set 6=1) and 

_ (4n(n - 1) + 12n)cn-2 - mgc^_i 

~ 4(n + 2)(n + 1) - 16 + 4(n + 2) ' ^ ^ 

with Co = 1 and Ci = —ml/20. At the horizon we use the variable y = x — 1 and the 
ansatz r^od/) = S^o^™?/" where we obtain from (p!6[) 

bn = 

(-4(n - l)(5n -2)+ ml- 16)6^-1 + (mg - 4(n - 2)(4n - 3))bn-2 - 4(n - 3)(n - 7)&„_3 

with bo fixed by the Taylor expansion at oo and bi = {m^ — 16)6o/8 and 62 = 
{ml — 48)6i/32 + mlbQ/32. We now find the mass eigenvalues by matching these 
two expansions. The results for ml < 1000 are shown in Table II below. 
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niQ 


niQ 


5.1102 


6.11 


7.8234 


9.35 


10.3591 


12.39 


12.8375 


15.35 


15.2909 


18.28 


17.7280 


21.20 


20.1528 


24.09 


22.5718 


26.98 


24.9868 


29.88 


27.3998 


32.76 


29.8088 


35.64 



Table II 



Values of mo for the O glueball in QCD^ obtained from matching 
the exact Taylor expansions at the horizon and oo in supergravity. rfiQ is 
the same mass normalized such that the lowest 0~^^ mass is 4.07. 



We observe that for the O three-dimensional glueball our exact values for the 
masses differ slightly from the ones obtained in the shooting technique in [^.0 As 
our first eigenvalue essentially agrees with the one in that reference we confirm the 
agreement between the ratio of the lowest mass O and O'^^ glueballs in QCD^ in 
supergravity and on the lattice reported in p8[. For completeness, in Fig.3 below we 
plot the wave function rjo corresponding to the mass eigenvalue mo = 27.3998. 



^However, the mass ratios are still in excelent agreement. 
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X 10" 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Fig3: This plot shows the wave function riQ for the O gluebah in QCD^ 
with mo = 27.3998. The plot is obtained from the Taylor expansions at the 
horizon and oo which agree perfectly. 



Finally, we will examine the 0~^^ glueball in four- dimensional QCD. The ap- 



propriate wave equation is in this case p8 



(19) 



with r = . At oo the Taylor expansion is of the form, with x = 1/r, riQ{x) = 
J2'^=o Cn.x'^'^^ where cq = 1, ci = — mg/lS and C2 = m^/792. The recursion relation is 

(4(n - - I) + 2{n - i))c„_4 - mgc^ ^^^^ 



C„-l 



4(n+|)(n + f)-10(n+|) 
At the horizon we use the variable y = x — 1 and the ansatz ?7o(y) = Y^'^=o bnV^ which 
gives 



bn 



+1 
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-2{n - 2){2n - 5)6„_2 - 2(n - 1)(3 + 8{n - 2))6„_i - (24n(n - 1) - + 6n)5„ 

12n(n + 1) + 12(n + 1) 

where 6o is fixed by matching with the Taylor expansion from oo and bi = —mlbo/12 
and 62 = ^-0(^0^6)60/576. The masses we obtained are exhibited in Table III below. 



/n++ ^2 




26 Q4Q8 


5 1 91 S 


63 8820 


7 9926 


114.1326 


10.6833 


177.7429 


13.3320 


254.7283 


15.9602 


345.0944 


18.5767 


448.8437 


21.1859 


565.9776 


23.7903 


696.4967 


26.3912 


840.4013 


28.9897 


997.6925 


31.5863 



Table III 

Values of ml and m^ for the O^^ glueball in QCD4 obtained from 
matching the exact Taylor expansions at the horizon and 00 in supergrav- 
ity. 



In this case our exact mass eigenvalues are also in close agreement with the ones 
presented in . 
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4 Conclusions 



We have examined the eigenvalue problems which feature in Witten's generaliza- 
tion of the conjecture by Maldacena regarding large supersymmetric gauge theory 
at high temperature and strong coupling. We have studied the eigenvalue problem 
through exact analytical expansions (at both the horizon and infinity) and evalu- 
ated these exact expressions numerically. We have analyzed carefully the behaviour 
of eigenfunctions at the horizon and discussed the boundary conditions. We have 
emphasized the fact that the correct criteria for selecting the wave eigenfunctions 
are normalizability and smoothness at the horizon, have shown that no smooth solu- 
tions exist with vanishing derivative (Neumann boundary condition) at the horizon 
and have given a construction of such smooth solution. Given that we are using ex- 
act analytic expressions for the various wave functions, our mass eigenvalues can be 
determined to any desired precison. 

Our values for the glueball masses are in agreement with the ones found in 
p8| and we have explained why the two techniques give identical results to this level 
of accuracy. Since we used exact formulas for the analytic expansions of the wave 
function solutions we believe that our work reinforces the good agreement between 
various glueball spectra obtained in supergravity and on the lattice as was already 
described in p8[. We hope that the results of this analysis will be of use for further 



studies of the conjecture. We expect that indeed they will be necessary as soon as 
higher precision in the mass values becomes important. 
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